On the use of Purcell factors for plasmon antennas 
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The Purcell factor is the standard figure of merit for spontaneous emission enhancement in microcavities, that 
has also been proposed to describe emission enhancements for plasmonic resonances. A comparison is made of 
quality factor, mode volume and Purcell factor for single and coupled plasmon spheres to exact calculations of 
emission rates. The paper explains why the Purcell factor is not appropriate for plasmon antennas. 
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It is a paradigm of quantum optics that the internal dy- 
namics of a quantum system can be controlled by placing the 
quantum system in a photonic environment that is resonant 
with radiative transitions of the sourest. Indeed, according to 
Fermi's Golden Rule, the spontaneous emission decay rate is 
proportional to the density of photon states that the photonic 
environment offers for spontaneous decay ^ A resonant mi- 
crocavity with a high quality factor Q, and a mode confined in 
an ultrasmall volume V enhances the density of photon states 
by a factor know as the Purcell factor Q/V^^^. Creating mi- 
crocavities with the highest Q/V is a central target in solid 
state quantum optics, with success being reported with pho- 
tonic crystals, micropillars, and Si02 microspheres and mi- 
crotoroids^'^. 

Recent advances in nanofabrication have provided new im- 
petus to the field of plasmonics, that creates highly confined 
electromagnetic resonances by using the response of free elec- 
trons in metals'*. In particular, there are widespread activities 
that exploit plasmon antennas for enhanced fluorescence and 
quantum optics with single plasmonsi"— . Due to absorption 
and radiative loss, the quality factor of plasmonic devices is 
extremely low {Q < 10^). Nonetheless, the truly nanoscale 
confinement leads to estimates for the Purcell factor of order 
F = lO'^, competitive with dielectric cavities""'-'. Motivated 
by the tremendous success of the Purcell factor as a quantita- 
tive figure of merit for dielectric cavities, several workers have 
also adopted this figure of merit for plasmonics'^"'*. Here, I 
demonstrate that the Purcell factor does not provide a quanti- 
tative gauge for spontaneous emission enhancements in plas- 
monic antennas. 

The spontaneous emission rate F of a dipole emitter (tran- 
sition frequency ui, transition dipole moment d) is governed 
by the local density of photonic states (LDOS) N{u!,r,ed) 
that enters Fermi's Golden Rule {cf. p. 273-275 ir^i) through 
r(w, r, ed) — V^^^ N{r, ui, e^). Here, r is the source position 
and e^ the dipole orientation. In terms of the electric field 
Green dyadic (G(r, r',uj) the LDOS equals^^ 



iV(a;,r,ed) 



6uj 



(ej •Im((G(r, r,w)) • e^) 



(1) 



This formulation is rigorous even for dissipative plasmonic, 
systems, in which case Eq. ([T]l includes quenching'**. The 
Purcell factor for dielectric cavities rests on two assumptions. 
First, in non-dissipative systems, normal modes can be de- 
fined, and the LDOS (dJ can be rewritten as a sum over all the 



normal modes 

iV(wo,r,ed) = 3^ |ej • Exir,u:)\^SiLJ - uo)- (2) 

The second assumption (see Gerard et al-') is that the sum 
over all modes in Eq. dU is dominated by just a single 
term. The (5-function is replaced by a Lorentzian spectrum 
of width Au = Lo/Q centered around the cavity resonance, 
and the spatial profile is set by the approximate eigenfunction 
of the cavity, normalized to contain a single photon within 
its mode volume V . Within this approximation, the LDOS 
-^cav('^cav, r,„ax,ed) at the spatial and spectral mode maxi- 
mum, normalized to the LDOS A'liost = naj^/(37r^c'^) of the 
host medium (index n) is well approximated by the Purcell 
factor 
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Maier has previously noted that reasonable numerical val- 
ues for F and V can only be obtained for plasmonic sys- 
tems if the energy density €\E\'^ that enters the mode vol- 
ume V = /e|i?pdr/max(e|i?p) is redefined'-' to avoid, 
e.g., negative energy densities'-'"'. The proper energy den- 
sity that is positive and real even when Re(e) < 0, and that 
accounts for energy stored in the material was derived by Rup- 
pjj^Bje pollowing Ruppins formalism."', the energy density 
used here is (Ree + 2ajlme/7)|i?p, appropriate for a Drude 
model e = 1 — ^^/^(a; + 17) for the dielectric constant of 
silver (ujp = 7.9 eV, damping rate 7 = 0.06 eV). 

In this work, the commonly employed procedure to obtain 
F described above is compared with exact calculations of the 
LDOS for the textbook case of metal spheres. First, I con- 
sider calculation of the mode volume. For dielectric cavities, 
one typically runs a numerical solver to obtain the local field. 
As long as the excitation has spatial and spectral overlap with 
the cavity mode, a high F ensures that the calculated spatial 
distribution of the local field of the resonance is almost inde- 
pendent of the exact excitation. The field distribution (total 
field minus excitation) is used to evaluate the mode volume 
integral, with integration limits truncated at the computation 
domain. Also for plasmon antennas, far field excitation yields 
near-field patterns with large local field enhancements, repre- 
sentative of the resonance mode profile. In this work, the exact 
Green function for a sphere'^ is used to calculate the local en- 
ergy density using plane wave excitation. FigurelTfa) shows a 
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FIG. 1 : (Color online) Contour plots of electric field \ES^ and energy 
density W for a 35 nm radius silver sphere in vacuum (plane wave 
excitation from below, horizontal polarization). The colorbar ranges 
from to 107 in units of the incident \E\ , resp. to 1 18 in units of 
W of the incident wave, (b) Mode volume J Wdr/raeix{e{r)\Ef) 
evaluated over a spherical integration domain truncated at R, plotted 
versus R. Vertical line: particle radius. Dashed line: linear diver- 
gence due to radiation loss. Horizontal line: mode volume, (c) Mode 
volume (red), polarizability (blue) and physical volume (dashed) ver- 
sus sphere radius for Ag spheres in vacuum. The right hand scale 
refers to the quality factor Q (green), (d) Red curve: Purcell en- 
hancement calculated from Q and V. Dashed: exact radiative decay 
rate. Solid curves: total decay rate for a dipole at 1 nm and 0.1 nm 
from the sphere. 



Strongly enhanced and confined field (enhancement relative to 
incident |i?p by ^ 10^) just outside the sphere (sphere radius 
35 nm). The energy density reaches a maximumjust inside the 
sphere. In FigurelTJb) a plot is shown of the mode volume in- 
tegral evaluated over a truncated integration domain of radius 
R around the plasmon sphere. Figure [TJb) exemplifies a hith- 
erto unnoted problem that is generic for all radiatively damped 
resonances'^, including all plasmon antennas. The mode vol- 
ume integral shows a linear divergence with the truncation ra- 
dius R. The divergence is understood by noting that radiative 
loss signifies constant integrated far-field flux, and hence a 
l/R^ asymptotic fall-off of |i?p. For dielectric cavities, this 
divergence is generally not apparent as high Q's imply a very 
small slope of the linear tail in Fig. [TJb), that is often negli- 
gible in a finite numerical computation domain. Leaving the 
fundamental treatment aside '^, I note that for dielectric cavi- 
ties a rigorous mode volume can be retrieved by subtracting 
the radiative part (linear divergence). Figure [TJc) shows the 
divergence-corrected mode volume for silver spheres of dif- 
ferent sizes. The mode volume is generally comparable to, or 
even smaller (large spheres) than the physical volume, since 
the field energy is mainly stored near the metal surface. 
To obtain the quality factor Q, Mie extinction cross sec- 
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FIG. 2: Radiative decay rate (dashed) and total decay rate (joined 
symbols) enhancement along the centerline joining two 25 nm silver 
spheres with a 10 nm gap. Enhancement is relative to the rate in the 
n = 1.5 host medium. Red: prediction from Purcell theory. The 
dipole orientation is along the centerline. Right: contour plot of |i?p 
(scale from to 1890 in units of incident \E\'^, dimer illuminated by 
a vertically polarized plane wave). 



tions'^ are analyzed. By dividing out the non-resonant cj^ 
scaling, Lorentzian line shapes are retrieved from which Q 
can be extracted. Figure [Tfc) shows that the quality factor de- 
creases from Q ^ 40 (limited by absorption) in very small 
particles to Q ~ 3 for larger particles. Combining Q and V in 
Fig- Ed) reveals predicted Purcell factors that range from 10^ 
for the smallest spheres to below 10^ for the 50 nm spheres. 
The key question is in how far this Purcell factor F correctly 
predicts the emission rate enhancement calculated rigorously 
from ImG'^. For this comparison, the maximum total decay 
rate enhancement almost at the particle surface (separations 
shown: 1 nm and 0.1 nm), where the mode field is maxi- 
mum, is plotted in Fig.fTfd). The Purcell factor clearly under- 
estimates the total radiative decay rate by an order of magni- 
tude. Strictly speaking the, Purcell factor is an estimate for the 
LDOS, i.e. the total decay rate. The Purcell factor estimates 
the radiative decay reasonably for particles > 15 nm radius. 

One might argue that single plasmon spheres set a poor 
example, as they do not provide very high local field en- 
hancements typical for plasmonics. Figure |2| presents an 
analysis of the Purcell factor for two spheres with a narrow 
gap, with a highly confined gap mode. The calculation uses 
a multipole expansion multiple scattering code provided by 
Garcia de Abajc"^ to analyze a dimer of two 25 nm radius 
silver spheres separated by a 10 nm gap, in a n == 1.5 host 
medium. A scattering resonance with Q ~ 5.7 occurs near 
A = 506 nm, with a 1890 times \E\'^ enhancement polarized 
along the dimer axis. Purcell analysis predicts a mode volume 
V = 1.32 • 10^^ iJ.nr' and a rate enhancement ranging from 
F r-^ 1300 at a sphere surface to about 700 times exactly in 
between the spheres. The rigorous calculation predicts a rate 
enhancement of ^ 1200 times exactly in the middle of the 
gap. Again the Purcell factor underestimates the actual total 
and radiative decay rate enhancements by almost a factor of 2. 

Referring to the two assumptions underlying the Purcell 
factor, in general both assumptions will break down: Firstly, 
a set of normal modes can generally not be defined for disper- 
sive absorbing media. Secondly, if an expansion in modes can 



be generalized (as in the case of multipole expansion for the 
Green function of a plasmon sphere) no single term is dom- 
inant, even if the resonance is clearly due to one term. For 
instance, at the dipolar plasmon sphere resonance, the LDOS 
is dominated by non-resonant coupling of the emission dipole 
to higher order multipole moments of the plasmon sphere**-^". 
This nonresonant coupling is not contained in a Purcell fac- 
tor estimate based on Q and V, and corresponds to immediate 
quenching of emission when very close to the metal. The total 
decay rate diverges as z^^, where z is the separation between 
dipole and metal surfac o^^'^" . Such quenching occurs gen- 
erally in both planar and curved plasmonic structures^", and 
unfortunately always exactly coincides with the field maxi- 
mum where one seeks to apply the Purcell factor. In the single 
sphere case, the fact that the Purcell factor reasonably approx- 
imates the radiative rate is consistent with the fact that most 
non-resonant coupling is quenching. For the dimer case, if the 
dipole is in the middle of the gap, multipole contributions en- 
hance both the total and radiative decay rate while maintaining 
a high quantum efficiency. 

Based on the above, one is led to conclude that the Purcell 
factor is not adequate for plasmon antennas. One might ex- 
pect that moving the emitter outside the regime of the extreme 
near field limit allows one to retrieve the Purcell enhancement 
prediction, taking into account that displacement and detun- 
ing of the emitter causes the enhancement to decrease from 



F in proportion to the mode field and the Lorentzian cavity 
lineshape''^. Even then there is a fundamental difference that 
sets plasmon antennas aside. Contrary to Purcell theory, the 
emission rate of a dipole near a nano-antenna can follow a 
dispersive Fano-type frequency dependence^'^^. As first dis- 
cussed by Dulkeith et al. — the lineshape can be understood 
as destructive interference of radiation emitted by the emis- 
sion dipole with that radiated by the plasmon antenna. This 
phenomenon occurs for all antennas where coupling to radi- 
ation is important. Incidentally, in this regime the decay en- 
hancement can be well understood in a point dipole model —. 
The only relevant volume that quantifies the interaction in this 
case is the polarizability (Fig. [Til, not the mode volume. 

In closing, the Purcell factor appears to be unsuited as a 
general gauge to classify the emission enhancements of plas- 
monic resonances quantitatively. While predicted Purcell fac- 
tors may act as a quahtative design guideline, in practice one 
always needs to confirm by a full calculation that non-resonant 
channels are negligible. 
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